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ABSTRACT 

We use a principal components analysis of radio-selected (3CRR, 6CE and 7CRS) AGN 
datasets to define two parameters related to low-frequency (151 MHz) radio luminosity 
L151 and [OIII] luminosity Lroiii] : a parameter a encoding the £151 — £[oiii] correlation 
and a parameter j3 encoding scatter about this correlation. We describe methods 
for constructing generalized luminosity functions (GLFs) based on a, f3, redshift and 
schemes for unifying quasars and radio galaxies. These GLFs ca n be used to generat e 
radio luminosity functions (RLFs) which improve on those of lWillott et"aIT<|200H) . 
mostly because they incorporate scatter and are therefore much smoother. 

Luminosity-dependent unified schemes (e.g. a receding-torus scheme) have been 
invoked to explain the low quasar to radio galaxy fraction at low a and the differences 
in emission-line luminosities of radio quasars and radio galaxies. With the constraints 
of the 3CRR, 6CE and 7CRS datasets and radio source counts, our GLF approach 
was used to determine whether a receding-torus-like scheme is required if there are 
two populations of radio sources: one at low a, consisting of 'starved AGN'; the other 
at high a consisting of 'Eddington-tuned AGN'. Because of the overlap between these 
two populations and the effects of the (3 parameter, schemes with or without a receding 
torus can produce a low quasar fraction at low a and differences in [OIII] luminosity 
between radio galaxies and quasars. The receding torus may be a physical process 
important in one or more populations of radio sources, but this is not yet proved 
either by the quasar fraction or the emission-line properties of radio-selected samples. 
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1 INTRODUCTION 

It has been recognised for some time that there is a strong 
positive correlation between the extended-radio luminosities 



and narrow-emission-line 


luminosities of 3C radio sources 


( Baum & Hcckman 


1989; Rawlings et al. 


119891: iMcCarthvl 



1993), and by combi ning the 3C sample with the 7C Red- 
shift Survey (7CRS), Iwillott et al.1 <1999ft found that these 
correlations were not primarily due to redshift. This sug- 
gests that the sources of narrow-line emission and radio jets 
are linked, with expla nations for this link rangi ng from the 
effects of e nvironment l|Du nlop & PeacoclJl993j), to jet-disk 
symb iosis jRawlings fc Saunders! Il99ll : iFalcke fc Biermannl 

Unified schemes for radio galaxies and radio quasars 
propose that they are the same objects viewed at differ- 
ent angles between their radio axis and their line of sight. 



An obscuring torus is invoked to hide the nucleus at large 
angles to the jet axis i.e. in radio galaxies. Emission from 
narrow-line regions is believed to be broadly independent 
of the jet-axis orientation, as it is emitted from a region 
larger than the extent of the torus. The opening angle of 
the torus 0trans marks the transition from the object being 
viewed as a quasar to a radio galaxy. We therefore expect 
for simple, arguably naive, unified schemes, where Otrans 
is constant, that the distribution of emission-line luminosi- 
ties should be simi l ar for radio quasars and radio galaxies. 
IJackson fc Browne! dl990^> found that [OIII] luminosities of 
radio quasars are ~ 5 — 10 times more luminous than ra- 
dio galaxies. I n similar comparisons of radio quasars and 
radio galaxies. iHes. Barthel_fc_F_osbjirj j[l993jl found no dif- 
ference using [QUI. and lJackson fc Rawlingsl il997t) found 
no difference in [OIII] at high redshift. In summary, some 
studies comparing the narrow emission line strengths of ra- 
dio quasars and radio galaxies are seemingly in agreement 
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with the predictions of the simplest unified schemes, whilst 
others are not. 

These seemingly contradictory, results can perhaps be 
understood in the context of a 'receding torus' model for 
AGN . This model (lLawrencelll99ll ; iHill. Goodrich fc DePovl 
1996) proposes that the inner radius r of the obscuring torus 
is determined by the radius at which dust sublimates, scaling 
as r oc Lpjfof Assuming that the half-height h of the torus 
is independent of L p hot, then the half-opening-angle of the 
torus, = tan -1 r/h, will be larger in the higher-L p h ot ob- 
jects. This means that more luminous objects, with higher 
Lphot , are more likely to be viewed as quasars so that quasar 
fraction rises systematically with Lphot • It also means that 
orientation-independent quantities which scale with L p hot 
will be higher on av erage for quasars than for radio galax- 
ies. lLawrencel dl99lT) found that the fraction of narrow-line 
objects in the low-frequency-selected 3CR sample decreased 
with increasing radio luminosity, and that narrow-line ob- 
jects have weaker [OIII] at a fixed radio luminosity. He ar- 
gued that this was inconsistent with the simplest unified 
schemes but that it could be explained by a cone- angle de- 
pende nce on luminosity, e.g a receding-torus model. ISimpsonl 
( 1998) argued that the [OIII] emission line is a much bet- 
ter indicator of Lphot than [Oil] , explaining why differences 
between quasars and radio galax ies are more obvious in 
[O III] than f OII]. ISimpsonl i2003f) revisited the arguments 



of ISimpson] (|l99£ ) and corrected a s mall error (com pare 
Fig. 2 of lSimpsonll2003l with Fig. 4 of ISimpson| [i998). He 
concluded that, allowing for a spread in h, the receding torus 
model predicts that radio quasars should be a factor of a few 
brighter in L p hot than radio galaxies in samples exhibiting 
a w ide range of quasar fra ctions. 

IWillott et all <2002al) found that radio quasars have 
higher submillimetre luminosities by a factor of ~ 4 than 
radio galaxies of the same radio luminosity and redshift, 
a factor which cannot be reduced below ~ 2 by accounting 
for possible synchrotron contamination. This result supports 
the idea that the simplest unified schemes, where Otrans is 
constant, are not an adequate description of the relation- 
ship between radio galaxies and quasars: submillimetre emis- 
sion comes from cool dust grains in optically thin regions, 
and therefore radiates isotropically; moreover, such emission 
could not be obscured by a torus even if it were emitted 
from regions close to the nucleus. Submillimetre luminos- 
ity is therefore an orientation-independent quantity which 
might scale closely with L p hot either because the cool dust 
is heated directly by the quasar or because it is heated by a 
starburst whose luminosity scal es with Lphot- We con clude 
that the submillimetre study of IWillott et al . (2002a) is m 
quantitative ag reement with the receding torus model of 
ISimpsonl J2003) . 

Other arguments for a receding-torus-like model include 
the evidence that the fraction of lightly reddened z ~ 1 3C 
quasars decreases with increasing radio luminosity in agree- 
ment with the higher fraction of lines of sight expected to 
graze the torus at lower L D hot (at fixed h) in the reced- 
ing to rus model iHill et alJl996tlSimpson. Rawlings fc Lacvl 
1999). Also the near-infrared spectral indices of quasars from 
the 3CR sample are correlated with luminosity, whereas the 
optical spectral indices are uncorrelated with the quasar lu- 
minosity or orientation, so that the strength of the 'big red 
bump', relative to the ionizing continuum, appears to be less 



in the more l uminous objects JSimpson fc Rawlingj Eoool: 
ISimpsonlEooil . 

However, there is significant scatter in the relation- 
ship between radio and emission-line luminosities and a re- 
ceding torus is not the only way of explaining the differ- 
ences in the emission-line properties of radio quasars and 
radio galaxies or the luminosity-dependence of the quasar 
fraction (the fraction of objects that show broad emission 
lines). Two-population radio luminosity functions (RLFs) 
have been used to provide a best-fit to the 3CRR, 6CE and 
7CRS radio source r edshift surveys and radio source counts 
(Will ott et al] 120011. It is possible that a two-population 
model with a simple unified scheme in one population, com- 
bined with the effects of scatter, could mimic the effects 
of a receding torus in producing both emission-line differ- 
ences between radio galaxies and radio quasars and the 
gra dual rise in quasar fraction with emission-line luminos- 
ity. IWillott et al] (|2000T) found a drop in the quasar fraction 
of AGN at low luminosity, postulating that their results are 
consistent with either the emergence of a second population 
of low-luminosity radio sources, which lack a well-fed nu- 
cleus, or a gradual decrease in Otrans with decreasing radio 
luminosity. 

Traditionally, radio galaxies have been divided int o 
two populations, FRI and FRII, dFanaroff fc Rilevl ll974T) . 
based on radio structure but correlated with radio lumi- 
nosity such that above log 10 (Li5i/WHz _1 sr _1 ) ~ 25, ob- 
jects are FRII, whereas below that critic al L151, they are 
FRI. IChiaberge. Capetti fc Celottil (120021) found that the 
FRII population is inhomogeneous and that not all of them 
can be unified with quasars. They find a low-radiative- 
efficiency accretion, weak or absent broad-line emission and 
a lack of a significant nuclear absorbing structure for weak- 
jet, low-ionisation narrow-line galaxies. For broad-line ob- 
jects and obscured high-ionization narrow-line galaxies they 
see or infer intense ionizing emission, powerful jets and a 
torus-like absorber. In low-radio-frequency-selected samples, 
the fraction of objects with observed broad lines changes 
rapidly from ~ 0.4 above log 10 (Li5i/WHz _1 sr _1 ) ~ 26.5 
to ~ 0.1 bel o w thi s. The less luminous population from 
IWillott et all feOOlft is composed of FRIs and FRIIs with 
weak/absent emission lines and their more luminous popula- 
tion of strong-emission-line FRII radio galaxies and quasars. 
A two-population model like this is motivated by the fact 
that the presence or absence of a quasar nucleus (as shown 
by emission-line strength) seems likely to be connected 
to the properties of a cen tral engine of the radio source. 
iKaiser fc Alexander! il997tl . and others before them, have 
argued that radio structure is influenced by the environ- 
ment. It seems likely that a high jet power is necessary for a 
highly collimated non-dissipative jet, thus all high luminos- 
ity sources are FRII, but as lower jet powers are reached, a 
jet is more likely to disrupt resulting in a FRI source. The 
exact environmental density and homogeneity would deter- 
mine the radio luminosity threshol d of FRI/FRII wit hin the 
lower luminosity population. iLedlow fc Owenl il996tl found 
that the FRI/FRII division is proportional to the square 
of the optical luminosity of the host galaxy, which is plau- 
sibly related to such environmental effects, although other 
causes have been proposed. This sim ple picture is shown 
to be incomplete by the discovery by iBlundell fc Rawlings] 
( 2001) of an optically powerful quasar with FRI radio struc- 
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ture. It also seems that there is evidenc e that some FRI 
sourc es in 3C have hidden quasar nuclei llCao fc Rawlingsl 
2003), and o ne clear example of a 3CR FRI with broad 
lines (SCSSGjSuru^s^n^^lJ^^Oa) is known. In this p aper 
we follow iHine fc Longairl Jl97Sft . lLaing et all Il994h and 
l.lackson fc W all l|l999ft in assuming that the FRI/FRII di- 
vision is strongly influenced by environmental effects and is 
much less fundamental than a division based on accretion 
properties of the central engines. 

It is well known (e.g. lDunlop fc Pea cock 1990) that the 
evolution in comoving space density of radio galaxies de- 
pends on both L151 and z. In two-population models this can 
be expres sed as differential d e nsity evolution of thes e pop- 
ulations IWillott et all 1200 ll) . ISnellen fc Best) (1200111 have 
criticized this approach on the basis of the artefacts it pro- 
duces when the populations join. 

The primary motivation of this work was to investi- 
gate in a quantitative manner the differences between the 
emission-line properties of radio-loud quasars and radio 
galaxies alongside the luminosity dependence of the quasar 
fraction. In this paper we will quantify the three-dimensional 
distribution in [OIII] emission-line luminosity L[oiii], 151 
MHz luminosity L151 and redshift z using a new method ap- 
plied to complete samples. A principal components analysis 
(PCA) is used to find the axis which causes the most differ- 
entiation between the objects, and it is possible that it can 
identify parameters which are more physically meaningful 
than the attributes of each data point (radio and emission- 
line luminosity) used in the PCA. We carry out a PCA and 
define a generalized luminosity function (GLF) in this new 
parameter space. The GLF is then used to make proper com- 
parisons between receding-torus and two-population models, 
allowing for cosmic evolution of the populations. 

The paper is organized as follows. In Sec. [5] the data 
are described. In Sec. |H] the principal components analy- 
sis is described. In Sec. 0] simple one-population GLFs are 
described to enable a comparison between unified schemes 
with and without a receding torus. The effects of using two- 
population GLFs are presented in Sec. and in Sec. |S| wc 
compare our results with previous studies and discuss possi- 
ble meanings for the new parameters found by the PCA. We 
assume throughout that Ho = 70kms _1 Mpc -1 , J7m = 0.3 
and Qa = 0.7. 



2 DATA 

2.1 Classification of radio sources 

Most theories postulate the existence of a central accreting 
supermassive black hole in all radio sources, so any divi- 
sion between radio quasars and radio galaxies is bound to 
be problematic. Operationally, however, the division seems 
reasonably clean. Previous prescriptions for making this di - 
vision (e.g. I Jackson fc Rawlingsl Il997t IWillott et all Il99d) 
have defined radio quasars as objects whose integrated opti- 
cal luminosities are dominated by point sources (rather than 
the host galaxy). For most other radio sources, defined as ra- 
dio galaxies, there is no direct detection of a compact nuclear 
source or broad emission lines, altho ugh of course bot h could 
be present but obscure d. Following IHine fc Longairl (I1979T) 
and lLaing et alJ dl994T> . most authors tend to sub-divide 



radio galaxies into those with high-excitation narrow-line 
spectra (HEGs) and those with low-excitation narrow- line 
spectra (LEGs). We follow Ijackson fc Rawlingsl il997f) in 
defining LEGs as objects with (rest-frame) [OIII] equivalent 
widths of < 10 A, [Oil] /[OIII] ratios > 1, or both. 

Amongst HEGs, there are a number of cases where clas- 
sification as a radio quasar is arguably more natural than 
classification as a radio galaxy. This number has tended 
to increase as observational methods have improved, and 
become more varied and sophisticated. Such objects, of- 
ten called 'weak quasars' can be sub-divided into a num- 
ber of distinct categories: (i) objects with unobscured, 
broad-line optical nuclei which are ins ufficiently luminous 
to outsh ine the host gala xy (3C382, ICohen et al\ |1999|: 
3C39 0.3, iPerez et alJll988t 5C7.17, 5C7.118. IWillott et alJ 
1998); (ii) objects with lightly veiled (dust obscuration 
Ay ~ 1 — 5) broad-line nuclei seen vi a broad wings on 
the Ha lin e (3C33.1, 3C67, 3C 268.3, l Laing et al.1 19941 ; 
3C22, 3 C41.ISimpson et al.lfl999l 3 C109. ICohen et aljll99a 
3C303, lEracleous fc Haloernl Il99l or via broad Paschen 
lines (3C184.1, 3C219, 3C223. iHill et al.lll996ft . via broad 
wings on rest-frame UV lines (3C325, Grimes et al. in 

prep), or inferred fro m nuclear poi nt sources seen in 

the thermal IR f3C65 [Simpson et alJll999l : 3C79, 3C234, 
ISimpson. Ward fc Wall 2000) or, at high spa tial resolution, 
with the HST f3C455. lLehnert et al.lTl999Tl : and (iii) ob- 
jects with broad lines apparent only in optically-polarized 
light, pr esumably scattered from a nucleus with A y ~ 10 
(3C265. iDev fc Spinradlll99rJ; ISimpson et al.lll99gl : 3C226, 
Idi Serego Alighieri Cimatti fc Fosburvll994h . The thermal- 
IR observations "o7 1slmpson^t^lTi^99r show, nevertheless, 
that material of high (Ay 10) visual extinction is needed 
to hide any bright nuclei in m ost HEGs. It is highly plausible 
(c.f Cyg A. iQgle et al"1ll997ft that some of these highly ob- 
scured nuclei will eventually show scattered (i.e. polarized) 
broad lines in deep spectropolarimetry. 

For the LEGs there seems little doubt that any 
quasar nuclei are of comparatively low bolometric lumi- 
no sity, but whether or not they are naked (as argued 
bv IChiaberge. Capetti fc Celottil 119991 or buried (as pre- 
dicted by the receding torus model) is still open to de- 
bate. There is one clear example of a broad-line, optically- 
compact nucleus in a LEG (3C38 6, ISimpson et alJ Il996l : 
IChiaberge. Capetti fc Celottilll999l) . 

With this bewildering variety of ways in which the terms 
radio quasars and radio galaxies can become confused, it is 
clear that any division into two classes needs to be care- 
fully explained. Fig.Qis a cartoon of a classification scheme 
based on both the intrinsic (unobscured) bolometric lumi- 
nosity of nucleus and the visual extinction to the nucleus, 
ft is worth emphasising that how one goes about meaning- 
fully separating radio quasars and radio galaxies depends on 
the nature of the question being asked. Here, as discussed 
in detail in Sec. 1, we are principally interested in compar- 
ing a two-population-like model in which LEGs have naked, 
but very low luminosity, nuclei with a receding-torus-like 
model in which LboI controls Ay so that LEGs necessar- 
ily have high Ay. The natural division for this problem, 
illustrated in Fig. separates objects with any direct evi- 
dence for a quasar nucleus, veiled or not and weak or not, 
from those in which there is no such evidence. In most fig- 
ures we will plot weak quasars (WQs), which are signifi- 
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Figure 1. A cartoon showing an idealized way of classifying the 
optical counterparts of radio sources. The visual extinction (Ay) 
towards the compact optical nucleus is plotted against the (un- 
obscured) bolometric luminosity (Lboi) of this nucleus (in units 
of the luminosity L* at the break of the galaxy luminosity func- 
tion). The areas of this plane are labelled with the types of object 
they are likely to contain. The shaded region shows the location 
of objects defined as radio quasars (which have Mg < —23), 
with lighter shading illustrating the location of 'weak quasars' 
(WQ) - note that this includes both weak and/or veiled objects 
and, rather arbitrarily (see text), includes some quasar nuclei 
(like 3C265) in which broad lines are detect ed only after scat- 
tering and excludes some objects (like 3C321. ICohen et al.ll999h 
although scattered broad lines are seen. The lower luminosity 
limit for weak quasars was derived from the lowest luminosity 
7CRS weak quasar. The location of the LEGs in this plane is 
strongly dependent on whether a two-population-like scheme or a 
receding-torus-like scheme provides a better physical description 
of the radio source population. 



cantly (Av ?S 1) veiled and/or low luminosity, with a sepa- 
rate symbol. Whether or not one defines objects with broad 
lines only in scattered light is a matter of debate because 
broad lines could, in theory, be scattered from highly ob- 
scured nuclei. Considering such objects, we somewhat arbi- 
trarily classify 3C226 and 3C265 as weak quasars - in the 
latter case because a thermal-IR nucle us with Av ~ 10 is 
seen iSimpson. Rawlings fc LacvMl999ft - and 3C321 as a 
radio galaxy becau se although scattered broad Ha is seen 

Cohen et all 1 1999ft . transmitted broad Paschen a is not 

Hill et alJll996D . 



2.2 Complete samples 

A complete sample consists of every radio source in a cer- 
tain area of sky brighter than a specified flux-density limit 
at the specified selection frequency. Ideally, all the sources 
are identified optically with a radio galaxy or quasar, the 
redshifts and emission-line luminosities are determined for 
each source and the sample sky area and flux-density limits 
are known. Radio (151 MHz) and e mission line ([QUID data 
were u sed from the 3CRR sample o f|Laing^Rjlgv^^ Longairl 
Jl983ft . the 6CE sample ( iRawlings. Eales fc Lacvll200li a 
revision of the 6C sample of lEaleJl98Bft and the 7CRS sam- 
ple JWillott et al.ll2002bl) . 

The 3CRR sample is selected with Si78 ^ 10.9 Jy 




Figure 2. The 151 MHz luminosity L151 versus redshift z plane 
for the 3CRR, 6CE and 7CRS samples described in Sec. IO and 
Table The various symbols identify radio quasars and radio 
galaxies from different samples: 3CRR quasars (filled squares); 
3CRR radio galaxies (open squares); 6CE quasars (filled trian- 
gles); 6CE radio galaxies (open triangles); 7CRS quasars (filled 
circles) and 7CRS radio galaxies (open circles). Weak Quasars 
(see Sec. 12.11 are shown by filled stars in both the 3CRR and 
7CRS samples. The dotted line shows the flux-density lower limit 
for the 3CRR sample, the dashed lines show the upper and lower 
limits for 6CE and the solid line shows the lower limit for 7CRS. 
The dot-dashed lines mark the L151 limits of the 7CRS sample 
at z = 0.25,0.5,1.0,2.0,3.0, as used in Fig.|3 



(5151 ^ 12.4 Jy assuming a spectral index of 0.8), and with 
a sky area of 4.23 steradians. Three sources are excluded, 
(3C231 because the radio emission is due to a starburst 
and not an AGN, and the flat-spectrum quasars, 3C345 
and 3C454.3 as they are only in the sample because o f 
Doppler boosting), leaving 170 sources dWirlott et all200lft . 
Of these, 39 are classified as quasars, 20 as weak quasars, 
and 111 are radio galaxies, (see Table Q for a summary of 
all samples). The radio galaxies consist of 88 FRIIs and 23 
FRIs, and unless otherwise specified, the term radio galaxies 
will henceforth refer to both FRIIs and FRIs. 

The flux-density limits of the 6CE sample are 2.0 ^ 
5*151 ^ 3.93 Jy, and the sky area is 0.103 sr, and 58 sources 
are used (as one object 6C1036+3616 is occluded by a bright 
star and is excluded without bias), of which 49 are radio 
galaxies and 9 are radio quasars. 

The 7C-I and 7C-II subsamples of the 7CRS are used, 
consisting of 37 and 39 sources, with 0.0061 sr and 0.0069 
sr and S151 ^ 0.51 Jy and Sisi 0.49 Jy respectively. Two 
sources are excluded: 3C200 because it is already included 
in the 3C sample and 5C7.230 as it is a flat-spectrum object 
which is only in the sample because of Doppler boosting. Of 
these objects, 49 are radio galaxies, 23 radio quasars and 2 
weak quasars (5C7.17 and 5C7.118). 

Even coverage of the radio luminosity-redshift plane 
would provide the best constraint with which to derive lu- 
minosity functions. However from Fig. we can see that 
the complete samples described above have a very uneven 
distribution on the L151 — 2 plane, because the samples are 
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Sample 


Area 
sr 


Si5l lower limit 

Jy 


S151 upper limit 

Jy 


Q 


WQ 


FRII RG 


FRI RG 


Total 


3CRR 


4.23 


12.4 




39 


20 


88 


23 


170 


6CE 


0.103 


2.0 


3.93 


9 





19 





58 


7CRS(I and II) 


0.013 


0.5 




23 


2 


14 


5 


74 


Total 








71 


22 


181 


28 


302 



Table 1. Table of the complete samples used to constrain the GLFs. 



flux-density limited and quite small. However, the 3CRR, 
6CE and 7CRS samples do have very high spectroscopic 
completeness and there are reliable redshift estimates for 
the small fraction of objects without spectroscopic redshifts 
JWillott. Rawlines fc BlundeUll2001ah . 



2.3 Emission-line data 

The emission-line data for th e 6CE sample are given in 
R^wJjngs^Eale^^^LacvJ (l200ll) and for the 7CRS sample in 



Willott et all (l2002bl) .For the 3CRR survey, the emission- 



line d a ta are taken primarily from Ijackson fc Rawlines! 
19971) . iHirst. Jackson fc Rawlines! J20031), R. Laing and 



J. Wall (priv. comm. ) following fWillott et alJ dl99£ 

Since the sources in the samples described above have 
redshifts in the range < z < 4, it is not possible to measure 
the flux of the same emission line in each object with opti- 
cal spectroscopy and the availability of near-infrared spec- 
troscopy is limited. The [OIII] A5007 line is chosen as the 
second luminosity in the model because it is an excellent 
indicato r of the streng th of the underlying non-stellar con- 
tinuum jSimpsonll998l) . and it is probably pro duced at radii 
beyon d those obscured by the torus. Note that lWillott et alJ 
( 1999) used [Oil] because it was the line with the most mea- 
surements. 

In cases where no [OIII] line flux was available, other 
narrow- line fluxes were used to estimate it. Most commonly, 
the other line flux used was [Oil]. There are 57 objects 
in the 3CRR sample with both [OIII] and [Oil] measure- 
ments. This enabled a best-fitting relation to be deter- 
mined of the form log 10 (L[ O n]) = a + & login (^[oiii]); where 
a = 4.74 ± 1.44, b = 0.86 ± 0.04 (see Fig. El . The 'best-fit' 
line was calculated using an algorithm that minimises the 
sum of the squares of the perpendicular distances from the 
data points to a line with an adjustable slope and intercept. 
One reason why this relation is not a proportionality may 
be the systematic changes of [Oil] /[OIII] ratio with narrow- 
line luminosity expected if the effec tive ionization param - 
eter changes as a function of L p hot dSaunders et alJ|l989l) : 
HEGs, which dominate at high L p hot, have a systematically 
lower [OII]/[OIII] ratio than LEGs, which dominate at low 
Lphot- The most common other lines used were Lyman-a, 
Mg II, [NelV], [NeV], CHI], Ha and CII] which were used 
to estimat e [OII I] using the average line ratios quoted by 
iMcCarthvl il993l) . 

Fourteen radio galaxies in the 3CRR sample and three 
radio galaxies in the 6CE had no emission-line data and 
these objects are identified in Fig. |3] These objects are the 
FRI radio galaxies 3C83.1B, 3C288, 3C296, 3C310, 3C314.1, 
3C315, NGC6109, NGC6251, NGC7385 and 3C465, the 
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Figure 3. [Oil] emission-line luminosity against [OIII] emission- 
line luminosity for those objects with data in both lines from the 
3CRR sample. The solid line shows the best-fitting relation. 



FRII LEGs 3C427.1, 6C1143+3703 and 6C1159+3651, and 
the FRII HEGs 3C68.2, 4C74.16, 3C292 and 6C1129+3710. 
Emission-line data for these objects was obtained by gener- 
ating a random value with a Gaussian distribution function 
with the mean given by Eqn.|I]and with the observed scat- 
ter. 

Three 3CRR quasars (3C9, 3C432 and 3C454), six 
6CE quasars (6C0824+3535, 6C0913+3907, 6C1148+3842, 
6C1213+3504, 6C1220+3723 and 6C1255+3700) and twelve 
7CRS quasars (5C6.5, 5C6.33, 5C6.34, 5C6.95, 5C6.160, 
5C6.237, 5C6.279, 5C6.287, 7C0808+2854, 5C7.70, 5C7.87 
and 7C0825+2930) had [OIII] emission-line strengths esti- 
mated by assum ing a rest-frame equivalent width of 30A 
Jffirst et al.ll2003l) . and the optical frequency power-law in- 
dex a op t = -0.44 for A < 5000A a nd a op t = -2.45 for 
A > 5000A JVanden Berk et alllioOll) . The distribution of 
quasars with estimates of iroin] from equivalent widths is 
found to be compatible with the distribution of quasars with 
real measurements in each of the three samples and overall. 

The line streng ths for 6C0955+3844 a nd 6C1045+3513 
were taken from IWillott et ail (l2002al) . A new spec- 
trum, redshift, classification and £[011] has been found for 
3C325 (Grimes et al., in prep). An upper limit for the 
[OIII] emission-line strength for 3C386 was taken from 
ISimpson et alJ il996ft . All objects that had only an upper 
limit were treated as if the upper limit was the [OIII] lumi- 
nosity, but this is not expected to have a significant effect 
on any result. 

A complete updated list of measured 3CRR 
line luminosities and their references are available at 
http:/ /www-astro. physics. ox. ac.uk/~cjw/3crr/3crr. html A 
table of the 3CRR, 6CE and 7CRS emission-line strengths, 
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from either measurements or estimates, is available at 
://www-astro. physics. ox. ac.uk/~sr/grimes. html 

2.4 Radio source counts 

The number of known radio sources at 151 MHz is very 
much larger than the sources in the complete samples 
with redshifts described above and is a valuable extra con- 
straint. The source counts have b een deter mined from the 6C 
jHales. Baldwin fc Warnerll988l) and 7C (iMcGilchrist et alJ 
1990) surveys, which have much larger sky areas than the 
6CE and 7CRS redshift surveys. The 7C source counts go as 
faint as 0.1 Jy, and so provide a low-flux-density constraint 
on the GLF, even in the absence of redshift information. 









Eigenvalue 


1.86 


0.14 


Proportion 


0.93 


0.07 






e/3 


log 10 L151 normalized 


1 

s/2 


l 

s/2 


log 10 L[oiii] normalized 


1 

72 


1 

V2 



Table 2. The eigenvalues X a and A^ and eigenvectors e a and 
ep in terms of the normalized components of log 10 L151 and 
log 10 Lqiu for the principal components analysis. 



http 



3 PRINCIPAL COMPONENTS ANALYSIS 

PCA is a statistical technique that has been extensively 
used in astrophys i cs for the analysis of spectral data (e.g. 
iFolkes et a!]ll999l : iMadewick et al]l2002l) and for reducing 
the dimensionality of problems. It finds how a set of prop- 
erties of a sample of objects are inter-related, by identifying 
sets of parameters that always correlate and are most dis- 
criminatory between each source in the sample. The strength 
of each new parameter is measured by the amount it con- 
tributes to the variance of the sample. Sub-samples of data 
are identified from any clustering in the space defined by the 
new axes of the PCA. 

There h ave been many treatments of PCA in the lit - 
erature fe.g. iFrancis fc Wills! Il999l: iMurtaeh fc Hecklll987t) . 
PCA searches for the best-fitting set of axes to replace the 
initial set of axes corresponding to the attributes of the data, 
using the squared deviation of the points from the axes as 
the goodness-of-fit criterion and by enforcing orthogonality. 
This process gives rise to an eigenvalue equation. The matrix 
formed from all the eigenvectors rotates the original basis to 
a new set of orthogonal axes. 

Practically, it is necessary to standardize the variables 
under analysis, so that the components of the normalized 
attributes X = {xij} are given by 



(1) 



SjVn 



where the index i runs from 1 to n = 302, the total number of 
objects in the 3CRR, 6CE and 7CRS samples, and j — 1, 2, 
so that ri = log 10 L151 , r2 = log 10 L[oiii], for each attribute 
f is the mean and s the standard deviation. 



3.1 Results of the PCA 

The PCA was performed and the average value of log 10 I/151 
was found to be 26.73 with a spread of 1.05 and the aver- 
age value of log 10 L[om] was 35.46 with a spread of 1.10. 
This gave rise to the eigenvalues and eigenvectors given in 
Table |21 As can be seen, the first eigenvector contributes 93 
per cent of the total scatter. Also, the principal eigenvector 
a shows that the dominant relation between log 10 L151 and 
log 10 i[oni] is that they are positively correlated, (with an 
equal contribution from their normalized components be- 
cause there are only two components), and the secondary 



relation is an anti-correlation, much smaller in magnitude 
than the main correlation. 

We have now defined two new axes on which to define 
a GLF. The projections of the attributes on these new axes 
are given by a and f3, where 

and similarly 

From Fig. 0] we can see that the first transformation of 
the data (top to the centre plot) involves a normalization 
(Eqn. so that the axes are expanded. The second trans- 
formation of the data (centre to bottom plot) is a rotation 
(Eqns. Inland [3J of the data. There seems to be a cutoff in a- 
space at ~ —0.10, so that below this value there is only one 
(weak) quasar. Also, the radio quasars are biased to lower 
values of /3, because they have higher values of i[oni] f° r a 
given value of L151 . 

Notice that the objects from the 3CRR, 6CE and 7CRS 
samples can all be described by the same principal compo- 
nents. Since it can be seen that for any flux-density-limited 
sample, L151 and redshift are correlated (see Fig.|5Jl, it might 
be expected that different samples may have different princi- 
pal components. Tests have shown that these components do 
not vary significantly for various combinations of samples. 
The most obvious conclusion is that the correlation between 
emission-line and radio luminosi ties is much more im portant 
than any L151 — z correlations dWillott et aljll999l) . It can 
be shown from the normalized components that the rela- 
tionship between L[oiii] and L151 can be expressed as 



(4) 



/ £[Qiii] \ _ „ , , 7 / Liu. 

For comparison, ISerieant et al.l <ll99ct) found that there 
was a radio-optical correlation in steep-spec t rum quasars 
of the form L phot oc L^ " 1 . IWillott et all l|l99S|) found 
a less steep slope of L p hot oc Ljj* 1 for 24 radio loud 
quasars from 7CRS, and attributed the discrep ancy at least 
in part to the fact that Ser jeant et alJ ljl998ri did not im- 
pose a limiting absolu te magnitude in their definition of a 
quasar. IWillott! l|200lft found a slope of 1.00 ± 0.04 in the 
log 10 Lroii] — l°gio -^151 relation. All these previous studies 
adopted J1m = 1.0, = 0.0. The differences between these 
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Figure 4. Top: The log 10 £[om] — l°gio ^151 plane for the 3CRR, 
6CE and 7CRS radio galaxies and quasars. The symbols are as 
in Fig. [3] for all figures. Large circles surrounding open symbols 
identify radio galaxies that have estimates of £[om] from the 
log 10 i[om] — l°gio L151 correlation with the observed scatter. 
Larger circles surrounding filled symbols identify quasars that 
have i[oin] estimated by assuming a value of 30A for the rest- 
frame equivalent width of the [OIII] line. Centre: The normalized 



5io -^[0111]" 



normalized '. 



1 -L151 plane. The eigenvectors cor- 



responding to the first principal component (solid line) and second 
principal component (dotted line) are also shown. Directions of 
increasing a and /3 are marked. Bottom: The distribution of radio 
quasars and radio quasars in the a — j3 plane. 



slopes can probably largely be attributed to the subtle dif- 
ferences between quasars and radio galaxies. 

Values of a and for all objects in 
the samples are available from the table at 
http:/ /www-astro. physics. ox. ac.uk/~sr/grimes. html 



4 GENERALIZED LUMINOSITY FUNCTION 

In this section, we will define generalized luminosity func- 
tions (GLFs) assuming that all of the radio sources are 
drawn from the same population, for comparison with two- 
population GLFs, which will be presented in Sec.|3 A GLF 
p(a, /3, z) is modelled as 



p(a,(3, 
where 

/(*) 



polO 



-p(a-a break ) 



exp 



and 



g{(3) = exp 



2z\ 



(P~Po) 

2aj 



(•») 



(6) 



(T) 



and a is the first principal component, j3 is the second prin- 
cipal component, op is the scatter, po is the normalization 
and z\ and zq define t he Gaussian evolutio n in redshift f(z) 
which fits the RLF in lWillott et alJ (l200ll) . The GLF is de- 
fined as the number of sources per unit a, per unit f3 per 
unit volume. Gaussian scatter in a direction perpendicular 
to the bes t -fittin g slope of log 10 L[on] — log 10 L151 is seen 
bv lWillottJ (1200 ll) . so Gaussian scatter in f3 is assumed. The 
function is a broken power-law in a, so that p — p\ when 
a > ftbroak and p = otherwise. The /3-offset parameter 
Po accounts for the fact that a GLF p(a,/3, z) derived from 
a symmetric distribution of objects in the (a, /3)-plane will 
produce a distribution that has an offset in the /3-direction 
for a flux-density limited sample. This arises because, for 
a given redshift, we select only the highest values of L151. 
From the definition of a, this means that we see only objects 
with low values of Ljoiii] for a given value of a, and thus 
we see objects with less negative values of /3 than would 
be seen in a sample that was not flux-density-limited. An 
analogous problem is that the distribution of sources with 
respect to spectral index in a complete radio-selected sam- 
ple varies with the frequency at which the sa mple is selected, 
ilKellermannlll964HWilliams fc Bridlelll967l) . 

We consider both a simple unified scheme with constant 
©trans and a receding-torus scheme with 



tan ©tra 



(£[oni]/£o) 2 tan9 , 



(8) 



where Lo is taken as the median value of Ltonii- 

A factor of sin 9 accounts for the probability of detecting 
an object inclined at an angle 6 to the line of sight. The 
detected object is a quasar if 6 < ©trans and a radio galaxy 
if 8 ^ ©trans- Therefore the GLF for radio galaxies is given 
by 



pKG(a,P,z) = p(a,f3,z) X COs6trans, 

and the GLF for radio quasars is given by 

Pd(a,P,z) = p(a,f3,z) X (1 - COs9trans). 



(9) 



(10) 
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We will denote the GLF with the receding torus as 1R and 
with the simple unified scheme as IS. 

A maximum-likelihood analysis was performed to opti- 
mize the free parameters of the GLF. For both GLFs, there 
are eight free parameters: ap,po,pi, «brcak ,/3 ,Zo,Zi and 
©o. In order to maximise the likelihood we must min imize 
S, wh ere S = — 2 ln(likelihood). Following iMarshall et alJ 



(11) 



1983), we find an expression for S, 



NQ 

! £ ] 

i=l 



-2 > hx[p{ai, Pi, Zi){l — cosBtr 



NRG 

»£ 

i=l 



+ 
+ 



ln[p(m,(3i, zi) cos 9tr ans ] 



WO— I J\ dzdlog 10 Li5idlog 10 L[oiii] 
dV 

PhgM— \ J\ dzdlog 10 Li5idlog 10 L[ O iii], 



where NQ is the number of radio quasars including weak 
quasars, NRG is the number of radio galaxies, Q — 
f2(log 10 L151 , z) is the sky area available from the samples 
for this value of z and Lasi, J is the Jacobian matrix 
for the change of variables from principal components to 
log 10 i[oni] and log 10 Li5i, and dz x dV/dz is the differ- 
ential co-moving volume element. In essence, there are two 
terms for quasars, and two for radio galaxies: one is the sum 
over all NQ quasars or NRG radio galaxies in the samples, 
and the last two terms are the integrals over the model being 
tested and should give ~ 2 NQ and ~ 2 NRG respectively 
for good fits. 

The errors associated with the parameters were found 
by calculating the components of the Hessian matrix (VVS*) 
at the location of the minimum, inverting this matrix to ob- 
tain the covariance matrix ([cr 2 ]ij = 2[(VVS') _1 ]ij), and tak- 
ing the lcr errors as given by the square root of the diagonal 
elements of this matrix (e.g. ISivialll996l) . 

To find relative probabilities for the GLFs withor with- 
out a receding-torus, the procedure set out bv ISivial dl996l) 
was used. The ratio of the posterior probabilities of GLF 1R 
with respect to GLF IS, Prs, is given approximately by 



Prs = 



P(lRjdata) 
P(lSjdata) 



5 lR|n 



Vdet(VVSis| 



x F, 



(12) 



N /det(VVS' 1R \S„ 



where det(VVSiR|s min ) is the determinant of the Hessian 
matrix for GLF 1R, evaluated at S — Smin. All of the free 
parameters are common to both GLFs, except for tan Oti-ans 
which is coded in a different way in each GLF (see Equation 
|5J but the prior ranges on this parameter are equal, giving 
F = 1. 



4.1 Results 

The maximum-likelihood routine was run for both the GLF 
with a receding torus 1R and the GLF with the sim- 
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Figure 5. The radio luminosity function at 151 MH z and 
z = 0.0001,0.25,0.5,1.0,2.0,3.0 from IWillott et alJ <200lft (solid 
lines) compared with the RLF from the optimized GLFs 1R 
(dashed lines) and IS (dotted lines). The circles show the location 
of the lower limit to the flux density of the 7CRS sample for each 
redshift, as indicated on Fig. |2] which shows the luminosity at 
each redshift below which the GLF is not well constrained. The 
z = 2 and 2 = 3 curves are identical, but the limit for z = 2 is 
above that for 2 = 3. 



pie constant-Qtra ns unified scheme IS. A downhill simplex 
( Pres s et al Il992ft method was used to locate the best-fitting 
parameters, and these are presented in Tabled] GLF 1R was 
found to be more likely than GLF IS by a factor of ~ 10 7 . 
The difference in likelihood between 1R and IS arises mainly 
from the differences in their unified scheme parameters, as 
most of the parameters describing the shape of the GLF 
are very similar. This naively implies that the simple non- 
luminosity-dependent unified scheme is strongly ruled out, 
as there exists an alternative physical model which is far 
more probable given the d ata. This has pu t on a quantita- 
tive basis the obvious (e.g. lLawrencelll99lf) statement that 
a constant-quasar-fraction model is inconsistent with the 
data. In Fig. 2] it is clear that the quasar fraction changes 

dra matically with 

IWillott et ail <200ll) derived three models for the RLF 
at 151 MHz from the low-frequency-selected 3CRR, 6CE 
and 7CRS samples. The GLF in a and /3 space, p(a, j3, z) 
can be easily transformed to give a RLF p(log 10 L151, z) for 
comparison 



p(^og 10 L 151 ,z) = 



p{a,p,z)\J\ ^dlo f 



10 L[oni]. 



(13) 



Fig. El shows the Model B RLF from IWillott et ail 

l)200lfl . which has been adapted to a VLm = 0.3, Qa = 
0.7, H = 70 km s" 1 Mpc" 1 cosmology, and the RLFs de- 
rived from the GLFs 1R and IS, which have almost exactly 
the sa me shape. Note that agreement with IWillott, et alJ 
(2001) is only expected to be good above points which rep- 
resent the flux-density limit of 7CRS. This is because only 
L15 1 — z plane data has b een used to constrain the GLF, but 
the lWillott et all (|200ltl RLF is also constrained by source 
counts. 

Simulations of the log 10 L[oiii] — l°gio ^151 pl ane from 
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log(po) 


Pi 


^brcak 


zo 


21 


lo gio(°"^) 


Po 


e 


1R 


-2.049 ±0.031 


27.08 ± 1.35 


-0.058 ±0.005 


1.951 ± 0.057 


0.550 ± 0.032 


-3.340 ± 0.048 


-0.021 ± 0.003 


42.0l 3 3 ; 9 9 


IS 


-2.060 ±0.097 


27.06 ± 1.73 


-0.058 ±0.006 


1.940 ±0.058 


0.545 ± 0.039 


-3.342 ± 0.051 


-0.020 ± 0.003 


45.8l 3 3 -° 



Table 3. Best-fit parameters for the Model 1R and Model IS. For the simple unified scheme, 0o is the constant transition angle, ©trans, 
whereas for the receding torus model, ©o is related to ©trans by Ean.lsl using the median log 10 (Lo) = 35.405 as a reference. The values 
of S m i n are 8701.89 for 1R and 8734.37 for IS. As the values of det(VVS) for the two models are comparable, 1R is more likely than IS 
by ~ 10 7 . 



GLFs 1R and IS are shown in Fig. and show that these 
simple models can reproduce the data reasonably well at 
high values of L151 and L[om] ■ O ne of the key features of the 
data, namely a lack of quasars at low values of L[oiii], is seen 
in the simulations from GLF 1R, whereas the simulations 
from IS clearly do not cause a difference in emission-line 
luminosities between radio quasars and radio galaxies. It is 
therefore easy to see why there is such a huge difference in 
likelihoods. 

We can compa re t he receding-toru s parameter with 
IWillott et ai] i200(J) and ISimpsonl (Il998f) . Fig. □ shows the 
quasar fraction as a function of L[oiii] for GLFs 1R and IS, 
as well as two-population GLFs which will be introduced 
in Sec.|5] It is clear that a constant-transition-angle unified 
scheme cannot reproduce the quasar fraction data from the 
3CRR, 6CE and 7CRS samples across the range of Lroiiii, 
while the one-population GLF with a receding torus pro- 
vides a reasonably good fi t to the data. Note th at the quasar 
fraction curve derived in lWillott et alJ i2000t) does not in- 
clude FRI objects. The same conclusions can be drawn from 
Fig. |H| where the quasar fraction is plotted as a function of 
a. The quasar fraction from GLF 1R also agrees reasonably 
well with the receding-torus type mod el parameters derive d 
from 3CRR and 7CRS L [ou] data from lWillott et al.N200rl . 

The scatter in the /3-direction, ap can be related to the 
scatter in the L151 — L[ iii] relation by the factor S2^2n, 
giving a value of 0.58, which is similar to t he scatter of ~ 
0.6 in the radio-optical correlation found bv lSerieant et alJ 
(1998), and t he scatt e r in t he L[oii] — £151 correlation of 
0.54 found bv lWillottJ feOOlft . 



5 A TWO-POPULATION MODEL 

The results of Sec. l4.1l implv that some mechanism is needed 
to cause the change in quasar fraction with a, recalling that 
a is a measure of the strongly correlated radio and [OIII] lu- 
minosities. An alternative mechanism to the receding-torus 
scheme could be a two-population model composed of a 
high-a population of radio galaxies and radio quasars re- 
lated by a unified scheme, in addition to a low-a population 
mostly composed of radio galaxies, but with a unified scheme 
that allows a small fraction of quasars, to account for ob- 
jects such as 3C386. This would certainly cause an increase 
in quasar fraction with a. 

A two-population GLF was created by adding together 
a low-a population GLF (pi) to a hig h-a population GLF 
(p 2 ). This method is similar to that of IWillott et al.1 (|200lTl 
for the construction of the 151 MHz RLF, and utilises the 
same evolutionary forms of both populations. Thus, 



GLF One-population schemes 



IS 
IP, 



Simple 

Receding torus 



Two-population Schemes 
GLF Population 1 Population 2 



2SS Simple 
2RS Simple 
2RR Receding torus 



Simple 

Receding torus 
Receding torus 



Table 4. Summary of the unified schemes used in each population 
for each GLF. 



Pl =p 10 10- pi(Q - Q - t ' l) exp [-10^-^.1)] g{P)h{z), (14) 
where g(/3) is defined in Eqn. J7J and 



fi(z) = {l + z) k for z< Zl 
= (1 + z 1 ) k for z > Z\. 



Also 



P2 = p 2 olO- p2(Q - Q - t - 2) exp [_io«( a -^)] g(/3)f2(z), (15) 



where 



Mz) = 



exp 
1 



--(- 

2 v 



for 
for 



Z < Z2a 



2^ z 2 b 

Z > Z2a- 

There is an exponential fall-off in the density of pop- 
ulation 1 for a > a cut and in the density of population 2 
objects when a < a cut . A factor q = s\s/n/\/2 is the fac- 
tor by which the range in a is smaller than the range in 
log 10 1/151, and acts to ma ke the exponential decline happen 
as in Will ott et all teOOll) . Fig.OOshows the two-population 
GLF. 

The first GLF we consider is a two-population scheme 
where the high-a population 2 consists of radio quasars and 
radio galaxies unified by a receding-torus scheme, and the 
low-a population 1 consists of radio galaxies and quasars 
with a simple unified scheme. Henceforth we will refer to this 
GLF as 2RS. We also consider the case where the both the 
low- and high-a populations are unified by a simple unified 
scheme and we will refer to this as the 2SS GLF. In addition, 
we consider a GLF for a two-population scheme which has 
quasars and radio galaxies unified by a receding torus model 
in both populations, which will be referred to as Model 2RR. 
The radio galaxy and quasar populations densities are given 
by 



PRG(a,/3, z) = pi COSe tl -an S1 ± p2 COS 6 t 



rans2 ■ 



(16) 
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Figure 7. The quasar fraction against [OIII] luminosity for GLFs 
1R (solid line), IS (horizontal solid line), 2RS (dashed line), 2SS 
(dotted line) and 2RR ( dot-dashed line). Als o shown is the re- 
ceding torus model from IWillott et alJ l200fl) converted to the 
correct cosmology and using [OIII] instead of [Oil], which was 
originally used (triple-dot-dashed line). The squares show the 
quasar fractions for the 3CRR, 6CE and 7CRS datasets in bins 



of log 



10 ^[Olll] 



luminosity. The error bars show the \/N errors. 




0.20 



Figure 8. The quasar fraction against a, for GLFs 1R (solid 
line), IS (horizontal solid line), 2RS (dashed line), 2SS (dotted 
line) and 2RR (dot-dashed line). The squares show the quasar 
fractions for the 3CRR, 6CE and 7CRS datasets in bins of a. 
The error bars show the vJV errors. 



Figure 6. The log 10 £[om] — l°gio ^151 plane for the 3CRR, 6CE 
and 7CRS radio galaxies and radio quasars (top) and simulations 
of log 10 L[oiu] ar >d logjQ I/151 data from the 1R (centre) and IS 
(bottom) GLFs. Symbols are as in Fig. |21 

PQ(a,(3,z) = pl(l - COSe tr an S1 ) + P2(l - COSe t rans 2 ). 

These GLFs are constrained by the 3CRR, 6CE and 
7CRS datasets as before (Sec. 14.111 and, because of the in- 
creased complexity of the two-population GLFs, we need the 
additional constraint of the 6C and 7 C source coun t s (de - 
scribed in Sec. E3J. In addition, the IWillott et alJ <|200ll) 
RLF had the constraint of the local RLF. Here we use 
the more recently determined local (z ~ 0.25) RLF for 



AGN from the 2dF Galaxy Redshift Survey dSadler et alJ 
2002). We simply use their normalization of the local RLF 
at log 10 (Li5i/WHz _1 sr" 1 ) = 24.0 (converted from 1.4 GHz 
to 151 MHz using a ra d = —0.8) as a prior. Without the use 
of this we would have very little constraint on the faint end 
of the GLF at low values of redshift. The slope p\ of the low- 
q population can simila rly be fixed at a value 6.4, derived 
using the results of ISadler et alJ i2002l) . 

The function to be minimised is now defined as 

Stotaj = Aln(xi) +Bln(xl) + Cln(xc), (17) 

where the parameters A, B and C are 'Lagrange mu ltipliers' 
or Bayesian 'hyper-parameters' dLahav et alJl200ol) . which 
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effectively control the relative weight applied to the data 
sets. The function to be minimised can be shown to be 



"total — N A ln(5) + iV B Hxl) + N c \n( X l RhF ) 



(18) 



where 5* is the likelihood function for two-population mod- 
els, analogous to 5* defined in Eqn. 11110 . Xsc is the \ value 
for the radio source counts and Xlrlf i s the \ 2 value for the 
local RLF. Na,Nb and JVq are the number of data points 
associated with the 3CRR, 6CE and 7CRS £[oiii] an d L151 
data (n = 302 objects), the binned source count data (30 
6C and 16 7C bins) and the local RLF constraint (1) respec- 
tively. The \ 2 ls given by 

i — nbins / \ 2 

2 \ ( fd&t&.i /"models 



E 



(19) 
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where dN = 2400(5 
counts. 



1.5 

min 



(20) 



Figure 9. The GLF, p(»,/3, z), integrated over (3 to produce 
p(a, z) for z = 1. The solid line is the GLF 1R and the dashed 
lines show the 2RR GLF. The total GLF and the contributions 
from population 1 and population 2 are shown. 



5.1 Results from the two-population GLFs 

The maximum-likelihood routine, where S is defined in a 
similar way to that in Eqn. 1111 was performed on the two- 
population models, to find best-fitting values of the GLF 
parameters, and the results are given in Table GLF 
2RR, with a receding-torus model in both high- and low- 
q populations was the best-fitting model but, in terms of 
probability, it is only preferred over 2RS by a factor of 
~ 2 (and over 2SS by a factor 1.3). The scatter about the 
log 10 i[oni] — l°gio relation was found to be 0.62. 

The form of the best-fitting GLF can be seen in 
Fig. |5] which shows p(a, z) (the GLF integrated over /3) 
against a. The two-population GLF is similar to the one- 
population GLF at high a, but they differ at low values 
of a, where the two-population model is additionally con- 
strained by radio source counts. The RLF derived from the 
two-populat ion GLF now agrees quite well with the Model 
B RLF from lWillott et al.l fcOOll) over t he full range of L 15 i 
and z (Fig. HUt . ISnellen fc Besd ll200ll) n ote that the two- 
population RLFs of IWillott et alJ (|2001j) have prominent 
humps due to the different evolutionary forms used for the 
high- and low-luminosity populations, which are identical 
in form to those used for the high-and low-a populations 
in this study. They argue that FRI and FRII radio galax- 
ies should not be treated as intrinsically different classes of 
object, and that evolution should simply be a function of 
radio power, so that the higher power FRII objects should 
undergo stronger evolution than FRIs. Since FRI objects 
are all drawn from the low-a population of the GLFs de- 
fined in this study, their evolution is in general weaker than 
FRIIs, which may be drawn from either population. The 
RLFs der ived from the GLFs are smoother than the RLFs 
found by IWillott et al.l i200ll) since at each value of L151 
objects are drawn from a band of objects in a. It seems that 
the RLFs derived from GLFs naturally produce smoother 
results because the effects of scatter, if only partially, arc 
taken into account. The 2RS, 2SS and 2RR two-population 
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Figure 10. The radi o luminosity function at 151 MHz, 
p(Li5i,z), from IWillott et alJ l200lf> (dashed lines) compared 
with the optimized two-population GLFs for model 2RS (solid 
lines), (2SS and 2RR give almost identical RLFs and are not plot- 
ted). The GLFs are plotted at z = 0.0001,0.25,0.5,1.0,2.0,3.0, 
and can be identified since the lines corresponding to each red- 
shift increase monotonically in normalization with redshift. The 
dotted line indicates the fixed slope adopted for the low-a pop- 
ulation, and the filled circle indicates the normalization of the 
local RLF at z ~ .25 and log^Lisi/WHz^sr- 1 ) = 24.0 from 
ISadler et alJ 120021) , Note the much less prominent 'hump' arte- 
facts in or new GLF-derived RLFs. 

models all give very similar GLFs, so that the small likeli- 
hood differences probably arise from the differences in their 
unified schemes. 

The best-fitting unified scheme parameters are com- 
pared with the quasar fraction as a function of £[0111] in 
Fig.|7|and as a function of a in Fig.|H| The effect of combin- 
ing a simple unified scheme with a two-population model can 
be seen, as the two-population GLF 2SS fits the data much 
better than the one-population GLF IS. The error bars on 
the quasar- fraction data are sufficiently large that, despite 
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^2b 


Po 
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Ooi 


e 02 


^2RS 


2RS 


-4.486 


-3.776 


31.47 


-0.0501 


0.0249 


0.437 


1.684 


0.447 


-0.0219 


-3.283 


3.436 


20.5 


29.9 


1.00 


2SS 


-4.486 


-3.776 


31.47 


-0.0500 


0.0249 


0.437 


1.684 


0.447 


-0.0219 


-3.283 


3.437 


20.6 


52.7 


1.59 


2RR 


-4.486 


-3.776 


31.47 


-0.0501 


0.0249 


0.437 


1.684 


0.447 


-0.0219 


-3.283 


3.436 


29.4 


31.0 


2.08 


2RR 
2SS 
2RR 


0.152 
0.143 
0.152 


0.020 
0.020 
0.020 


1.33 
1.34 
1.34 


0.023 
0.022 
0.024 


0.0015 
0.0015 
0.0015 


0.045 
0.035 
0.047 


0.044 
0.043 
0.044 


0.032 
0.032 
0.033 


0.0007 
0.0007 
0.0007 


0.019 
0.019 
0.019 


0.556 
0.516 
0.563 


21.7 
48.2 

+48.5 
-25.5 


+22.9 
-15.8 
17.3 
+24.4 
-16.5 





Table 5. Best-fit parameters for the two-population models. The upper section gives the optimized parameters and the lower section 
gives the errors on these best-fit parameters. 0oi and ©02 give the constant transition angles ©transx and 0trans 2 f° r a simple unified 
scheme in population 1 or 2 i.e. population 1 of 2RS and 2SS and population 2 of 2SS. For receding-torus schemes, ©oi arltl ©02 are 

related to ©transi and ©transn 

by Eqn. using the median log 10 (io) = 35.405. The values of S m i n , where S m i n is the minimum value 
of the likelihood function 5, are 2933.78,2933.91 and 2934.10 for 2RS, 2SS and 2RR respectively and det(VVS 2 Rs|s min ) = 2.12 X 10 43 , 
det(WS 2 Ssls min ) = 7.42 X 10 42 and det(WS 2 R R |s min ) = 3.58 X 10 42 . Inserting these values into Equation [T2l we find the probability 
relative to model 2RS, P 2 r,s- 
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Figure 11. The source counts from the 6C catalogue (stars) 
and 7C (squares), with the integrated source counts from the 
IWillott et ah! feOOll) RLF with contributions from the high and 
low populations (dashed lines), and from model 2RS (solid lines) 
with the contributions from population 1 and population 2 AGN 
shown separately. 



the very different shape of the quasar fraction curves, the 
non-luminosity-dependent schemes used in 2SS do not give 
a significantly worse fit than GLFs employing a receding- 
torus model in one or both populations. These GLFs, 2RS 
and 2RR, fit the quasar fraction data very well, differing 
only at low values of a, where 2RR predicts that the quasar 
fraction goes to zero at very low a. In comparison, 2RS and 
2SS seem to over-predict slightly the number of objects at 
low a, in order to be able to account for objects like 3C386. 
It is obvious in this context why there is little difference be- 
tween the likelihoods of the GLFs. 2RR can fit the second 
lowest data point in Fig.|7|and Fig.|H] whereas 2RS and 2SS 
can both fit the lowest data point. None of the GLFs can fit 
the third lowest quasar fraction data point in either Fig. [7| 
or Fig. |H] which seems to arise from the presence of weak 
quasars at approximately —0.1 < a < —0.3 (see Fig. 3J. 

The source counts derived from the GLFs are com- 
pared with the 6C and 7C source count data (Sec. 12.41 and 
the source counts derived from the RLF from lWillott et ail 



Figure 12. The redshift distribution dN/dz at a flux-density 
limit of 0.1 Jy for both po pulations of model 2 RS (solid lines) 
and of Model B RLF from IWillott et alj l200lj> (dashed lines). 
The plots are normalized such that the total number of sources 
in the sample is 1000. 



fcOOlf) in Fig. El The RLFs from IWillott et all (1200 ill do 

not give a very good fit to the source count data (re- 
duced x 2 ~ 3.17), under-producing the 6C counts at ~ 
0.7 — 3.0 Jy. However, the employment of the Bayesian 
'hyper-parameters' in the definition of the likelihood func- 
tion (Eqn. I18p leads to a sufficiently strong source count 
constraint, so that the source counts derived from our GLFs 
result in an improved fit to the data (reduced \ 2 ~ 1-68). 
The difference between the data and the RLFs can be ac- 
counted for by shot noise and by deviations of the number 
of objects in the 7CRS sample from the average number 
of s ources in an area of sky due to large-scale structure, 
(e.g. iBrand et alJl200fl . Another consequence of the tighter 
source count constraint is that, despite the good agreement 
of the RLFs at most redshifts (Fig. 1101 . the low-a population 
1 provides les s of the source count s than the low- luminosity 
population of lWillott et al.l i200ll) . and correspondingly the 
hig h-q population 2 pro vides more of the source counts than 
the lWillott et all (1200 ll) high-luminosity population, (which 
is also evident in Fig |12t , 
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These subtle changes can have significant impact on the 
redshift distributions predicted for new radio source redshift 
surveys. For example, the peak of the model redshift distri- 
bution is shifte d to lower redshifts for the GLFs compared 
to the lWillott et al.1 (1200 ll) RLFs at the flux density limit of 
a redshift survey at 0.1 Jy fFig 1121 . 

The simulations of the data for the two-population 
GLFs are shown in Fig. 1131 All three GLFs show that ra- 
dio quasars are more likely to be found at high values of 
L[oiii] compared to radio galaxies, although the distribu- 
tion of quasars in the high-a population is more evenly dis- 
tributed in 2SS because of the constant opening angle of 
the torus. The most obvious problem with the simulations 
is that all of the GLFs under-produce the number of objects 
at the lowest radio and emission line luminosities. This is 
a consequence of t he reduced normaliz ation of population 
1 compared to the Will ott et alJ (1200 if) RLF at 3CRR flux 
densities (see Fig. 1111 and Fig I12H . This could be fixed by 
invoking a more complicated model. 



6 DISCUSSION 

6.1 Quantitative Comparison of Unified Schemes 

A preliminary investigation (Sec.0J showed that a receding- 
torus scheme is strongly preferred over a simple non- 
luminosity-dependent unified scheme. However when a two- 
population model was investigated, there was little evidence 
to favour a receding torus in the high-a population. It is clear 
that, with enough scatter in the radio - optical relation, a 
small quasar fraction in the low-a population combined with 
even a constant value of #trans in the high-a population will 
mimic the rise of quasar fraction with emission-line lumi- 
nosity predicted by a luminosity-dependent unified scheme 
such as the receding-torus model. A simple unified scheme 
alone will not generate the observed increase in quasar frac- 
tion with a and the emission-line differences between radio 
quasars and radio galaxies but a receding-torus model or 
a two-population scheme with any sort of unified scheme 
in the high-a population will produce these effects. On the 
basis of quasar fractions and emission-line differences alone, 
there is not enough evidence to distinguish between possible 
luminosity dependences of the unified schemes, apart from 
ruling out schemes with no luminosity dependence, such as 
mo del IS. 

ISimpsonl (|2003h reviews the evidence for the receding- 
torus model highlighting the quasar fraction and emission- 
line difference arguments. However we have shown here 
quantitatively that the quasar fraction and the differences 
in emission-line luminosity in radio-selected samples are not 
very significant either individually or when considered to- 
gether. 

A major cause of the inability of the data to discrimi- 
nate between the competing models is that the quasar frac- 
tions of the 6CE sample is very different (15.5%) to the 
3CRR (34.7%) and 7CRS samples (33.8%). Fig. El shows 
the quasar fractions as a function of the limiting flux den - 
sity for the 3CRR, 6CE, 7CRS and 7CQ jRilev et alJll999T) 
samples and the curves predicted from the 2RS, 2SS and 
2RR models. This illustrates the fact that although the 2SS 
model leads to very different predictions than the receding- 
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Figure 14. The quasar fraction against the flux density at 151 
MHz Si5i from model 2SS (dotted), 2RS (dashed) and 2RR (dot- 
dashed). Also shown are the quasar fractions at the flux-density 
limits from the 3CRR, 6CE and 7CRS radio-selected complete 
samples used as constraints in the modellin g. The lower limit on 
the quasar fraction from the 7CQ sample iRilev et alJll99flh is 
also shown. The error bars are the \/N errors. 

torus model, they cannot yet be clearly distinguished by the 
data. 

To re-investigate the issue of the emission-line differ- 
ences between radi o quasars and radio galaxi es, we repeat 
the experiment of Ijackson fc Browne! lll990t) . They com- 
pared each of 12 3C radio galaxies with radio quasars with 
0.19 < z < 0.85, matched in redshift to within 0.05 and 
in luminosity to within 30 per cent, and find that quasar 
[OIII] luminosities exceed those of radio galaxies by a factor 
of 5 — 10. We plot histograms of the number of 3CRR ra- 
dio quasars and radio galaxies with 0.2 ^ z ^ 0.8 in L[oiii] 
bins, Fig. 1151 There is clearly an offset in the distributions, 
with the median value of £[om] being 36.10 for quasars and 
35.50 for radio galaxies, giving a shift of 0.60, so that the 
median quasar has 4 times more luminous [OIII] emission 
than the median radio galaxy. We found the distributions in 
Lyoiu] of 3CRR objects in this redshift range for the three 
two-population models. All three models reproduce a shift in 
the median value of £[om] between quasars and radio galax- 
ies. The 2RS model gave a shift of 0.78 compared to 0.56 
for 2SS and 0.64 for 2RR. This analysis shows that any of 
the two-popul ation models considered h ere can reproduce 
the results of [Jackson fc Browne! lll990t) so that their test 
does not provide clear support for the receding-torus model. 
All models, with or without a receding torus, predict that 
quasars are brighter than radio galaxies because the proba- 
bility of an object being classified as a quasar increases with 
a and L[onrj- 

6.2 The meaning of a and /3 

It is instructive to consider, if only briefly, possible in- 
terpretations for the axes found by the principal compo- 
nents analysis in the context of the receding-torus scheme 
and the two-population scheme. Since a accounts for the 
overwhelming majority of the scatter, it must be a man- 
ifestation of some important physical process. Higher val- 
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Figure 13. The [OIII] emission-line luminosity against the 151 MHz luminosity data and simulations of this data from the two-population 
GLFs for models 2RS, 2SS and 2RS. Symbols are as in Fig.0 



ues of a mean higher low-frequency radio luminosity and 
higher emission-line luminosity, i.e. more powerful objects. 
Emission-line luminosity is often assumed to be dominated 
by photo-ionization from a strong nucleus, though there is 
evidence that some part, specifically l ow-ionization lines like 
[Oil] , is powered by radiative shocks lllnskip et al 120021, . In 
this context an increase in accretion rate will cause an in- 
crease in both radio luminos ity and emission-line luminosity 
jRawlines fc Saunderil99lh . although it is not immediately 
obvious in what proportions. We tentatively suggest that a 
is an indicator of accretion rate and j3 is the scatter about 
this relation, which can probably be attributed to a number 
of factors: for example, environment and black hole mass. In 
this context and for the two-population model, it seems pos- 
sible that a C ut,i and a cu t,2 represent respectively a critical 
accretion rate above which AGN are shining at a significant 



fraction of their Eddington luminosity, and a critical rate 
below which they are shining at a tiny fraction. 

In essence, it can be shown that, for the very similar 
values of scatter in the radio and emission-line distributions 
as found for this sample, 

/ Li5iL[oni] \ , , 

a oc log 10 f , (21) 

\M51,av-fj[OIII],av / 

and 

/3oclog 10 f V (22) 

\M51,av/-k[OIII],av / 

iisi gives an indication of the jet power which is in 
turn powered by the central engine and £[om], arising from 
photo-ionization, is also a measure of the power of the cen- 
tral engine, a represents the product of the radio (mechani- 
cal) and optical (radiative) output, compared to the average 
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Figure 15. The top panel shows the number of 3CRR radio 
quasars (solid lines) and radio galaxies (dotted lines) at 0.2 §C 
z $C 0.8 at a given £[om] luminosity. The other panels show the 
distribution of 3CRR radio quasars (solid lines) and radio galaxies 
(dotted lines) at the same redshifts drawn from simulations of 
the 2RS, 2SS and 2RR GLFs. The vertical solid lines show the 
positions of the median value of I/[oin] f° r radio galaxies and 
radio quasars. 



product, (3 represents the ratio of radio output to optical 
output compared to the average ratio. 

As (3 probably represents the scatter from a number of 
different sources, it would be interesting to see if f3 correlates 
with any of the mos t likely contributions tow ards scatter. It 
has been shown bv iMagorrian et alJ lll99ct) that the black 
hole mass correlates with the mass of the bulge in galax- 
ies, and assuming a constant mass-to-light ratio, we expect 
the stellar luminosit y to correlate with the black hole mass. 
IWillott et all £2003) find that the best-fit K - z relation 
from the 3CRR, 6CE, 7CRS and 6C* samples is given by 



K = 17.37 + 4.53 log 10 z - 0.31(log 10 z) 



(23) 



In Fig. 1161 a measure of the stellar luminosity AK (the 
excess K luminosity over the K — z relation of Eqn. 1231 x ) 
is plotted against j3 for 160 radio galaxies from the 3CRR, 



1 IWillott et alJ 120031) show that the empirical K — z relation of 
Eon. l2lfl lies very close to the locus expected for a giant elliptical 
of fixed mass whose stellar population formed at high redshift and 
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Figure 16. The excess K magnitude s over the be st fitting K — z 
relation of Eqn.|23] AK, from lWillott et alJ <2003ft . against f3 for 
radio galaxies from the 3CRR, 6CE and 7CRS (I and II) samples, 
with symbols as in Fig. |2] Large circles surrounding two 3C ob- 
jects indicate objects known to be in clusters lHardcastle et all 



6CE and 7CRS (I and II) samples (IWillott et alJl2003l) . The 

Spearman rank correlation coefficient is —0.064 with signif- 
icance 0.42 implying that /3 is not significantly correlated 
with black hole mass. The correlation coefficient between a 
and K is —0.13 1 with significance 0.10. 

lLacv et alJ i200 ll) show that, for a sample of quasars 
from the FIRST Bright Quasar Survey (FBQS) and 
optically-selected quasars from the Palomar-Green survey, 
M BH ±0 ' 2 ( L / L Edd) 10 , where M BH is the black hole 
mass. This agrees with the slope of 1.05 found here in the 
logio ^151 - log 10 L 



10 -Mom] 



relation if the narrow emission-line 



luminosity, L[oni] is proportional to the bolometric lumi- 
nosity, L, and if the range in black-hole mass is small. We 
conclude that the scatter ap is not dominated by the black- 
hole mass. 

There are of course a variety of other factors which 
m ight correlate with 0. For ex ample, implicit in the work 
of lRawlings fc Saunders] (Il99ll) is that (3 might be expected 
to correlate positively with the richness of the radio source 
environment as, for a given jet power (assumed to scale 
with bolometric quasar luminosity), simple physical mod- 
els predict that radio sources have hig her radio luminosi- 
ties i n dense-gas environments (see also lBarthel fc Arnaudl 
1996), leading to more positive values of (3. Sadly, despite 
the existence of XMM-Newton and Chandra, there are as 
yet no comprehensive studies of radio source environments 
that can be used to make such a test. We label on Fig. 1161 
several radio sources know n to be in rich galaxy clusters 
from ROSAT observations jHardcastle et alJll98ct) . The ra- 
dio galaxy 3C295 (/3 = 0.054) is one of the highest-/? ra- 
dio galaxies in the sample, in agreement with the pre- 
diction, but 3C220.1 (which is not in Fig. 1161 as it does 
not have a measured K magnitude) has j3 = 0.015 and 
3C280 has f3 = —0.013, compared to the average value 



is evolving passively. AK represents (in magnitudes) the excess 
luminosity. 



16 J. A. Grimes, S.Rawlings & C.J.Willott 



of j3 = 0.002 5 for radio galaxies in the redshift range 
0.4 < z < 1. lHardcastle et alj l)l988ft also finds that five 
3CRR quasars with 0.35 < z < 0.75 reside in rich clus- 
ter environments, and they have an average f3 = —0.0024, 
compared to an average of —0.0037 for quasars in this red- 
shift range. So there are hints that this correlation may 
be present, but obviously the numbers involved are very 
small. We have also failed to find any significant correla- 
tion between j3 and projected radio source linear size D. 
Naively one might expect a negative correlation because 
synchrotron and adiabatic losses should cause a system- 
atic decline in radio luminosity over the lifetime of the 
radio source i f the .jet power and narrow-line luminosity 
remain fixed l|K aiser. Dennett- Thorpe fe AlexandeJ Il997t 
iBlundell. Rawlings fe Wirlottlll999D .' However, as discussed 
bv TwillTiT^^t^lTi^99Tr . such a correlation may be masked 
by a number of effects including changes with time of jet 
power and L p hot, and the boosting of emissio n-line fluxes 
in sm all sources due to radiative bow shocks jlnskip et alJ 




7 CONCLUSIONS 

A new approach to investigating unified schemes has been 
presented, based on a principal components analysis of the 
3CRR, 6CE and 7CRS complete samples. Generalized lumi- 
nosity functions have been derived based on these principal 
components: a which encodes the L151 — I/join] correlation 
and (3 which encodes the scatter about this correlation. The 
main advantage of this new approach has been that the uni- 
fied scheme parameters have been found by taking into ac- 
count the intrinsic scatter in the L151 — I/[oiii] correlation. 
The main conclusions to be drawn from this analysis are as 
follows. 

(i) A receding-torus model is strongly favoured over a 
simple non-luminosity-dependent unified scheme for GLFs 
with one population of radio sources. 

(ii) With the extra constraint of 6C and 7C radio 
source counts and the normalization of the local RLF, 
two-population GLFs were derived. The GLFs give rise to 
RLFs which are ve ry similar but smoother than those of 
IWillott et alJ i200lf) . There is very little difference in like- 
lihoods between a GLF with a receding torus in both pop- 
ulations, a GLF with a torus opening angle that does not 
vary with ionizing luminosity in both populations, and a 
GLF with a receding-torus in the high-a population and a 
constant-opening-angle torus in the low-a population. 

(iii) Two-population models reproduce the radio survey 
data well, and can provide a natural explanation for the 
rise in quasar fraction with emission-line luminosity and the 
emission-line differences between radio quasars and radio 
galaxies. The receding-torus may be a feature in both pop- 
ulations but this is not yet proved. 
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